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ABSTRACT 

We study the relation between the mass accretion rate, the jet power, and the black hole mass 
of blazars. To this aim, we make use of the Sloan Digital Sky Survey (SDSS) and the 1 1 
months catalog of blazars detected at energies larger than 100 MeV by the Large Area Tele- 
scope (LAT) onboard the Fermi satellite. This allows to construct a relatively large sample of 
blazars with information about the luminosity (or upper limits) of their emission lines used as 
a proxy for the strength of the disc luminosity and on the luminosity of the high energy emis- 
sion, used as a proxy for the jet power. We find a good correlation between the luminosity 
of the broad lines and the 7-ray luminosity as detected by Fermi, both using absolute values 
of the luminosities and normalising them to the Eddington value. The data we have analyzed 
confirm that the division of blazars into BL Lacs and Flat Spectrum Radio Quasars (FSRQs) 
is controlled by the line luminosity in Eddington units. For small values of this ratio the object 
is a BL Lac, while it is a FSRQs for large values. The transition appears smooth, but a much 
larger number of objects is needed to confirm this point. 
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1 INTRODUCTION 

The classic division between Flat Spectrum Radio Quasars (FS- 
RQs) and BL Lac objects is mainly based on the Equivalent Width 
(EW) of the emission lines. Objects with rest frame EW> 5 A are 
classified as FSRQs (see e.g. Urry & Padovani 1995). Marcha et 
al. (1996) and Landt, Padovani & Giommi (2002) discussed the Ca 
H&K 4000-A (rest frame) break as a criterion helping to distin- 
guish BL Lac objects from low-luminosity radio galaxies. Further- 
more, Marcha et al. (1996) proposed that objects with a weak Ca 
break and with EW even larger than 5 A should be classified as 
BL Lacs. On the other hand, Scarpa & Falomo (1997) showed a 
continuity between BL Lac and FSRQs concerning the luminosity 
of the Mg II line, taken as an indication against a clear separation 
of blazars in the two subclasses. Landt et al. (2004), instead, con- 
sidered narrow lines, such as [O II] and [O III], and found that it 
is possible to separate intrinsically weak and strong line blazars in 
the [O II] and [O III] equivalent width plane. 

The classification scheme based on the EW of the broad lines 
has been adopted both because it is observationally simple, and 
because it was thought to measure the relative importance of the 
non-thermal jet emission over the thermal one. However, we now 
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know that the jet electromagnetic output is often dominated by the 
emission at higher energies (hard X-rays and 7-rays), and there- 
fore the EW of the optical emission lines is not a good measure 
of the jet dominance. Furthermore, the jet flux is much more vari- 
able than the underlying thermal emission, causing the measured 
EW to vary. Occasionally, a blazar with very luminous emission 
lines, that should be classified as a FSRQ, can instead appear as BL 
Lacs when the optical jet flux is particularly strong. Conversely, a 
BL Lac in a particular faint state could show broad emission lines 
that, albeit weak, can have EW greater than 5 A. In Ghisellini et al. 
(201 1 ; Gl 1 thereafter) we have therefore proposed a more physical 
distinction between the two classes of blazars, based on the lumi- 
nosity of the broad emission lines measured in Eddington units: 
LBh-R./L Edd . We proposed that when LBhn/L Edd >5 x 10" 4 the 
objects are FSRQs, and are BL Lacs below this value. Normalizing 
to the Eddington luminosity ensures the appropriate comparison 
among objects of different black hole masses. 

The sample of blazars studied in G 1 1 was limited, since it was 
based on a small sub-sample of bright FSRQs detected in 7-rays 
by the Fermi satellite during the first 3 months of operation (LBAS 
sample, Abdo et al. 2009), and on BL Lac objects detected by Fermi 
during the first 1 1 months (1LAC sample, Abdo et al. 2010a), with a 
relatively steep 7-ray energy spectral index a 7 (a 7 > 1.2). These 
BL Lacs occupy the region of the spectral index - 7-ray lumi- 
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nosity (a 7 -L 7 ) plane occupied mainly by FSRQs (see Ghisellini, 
Maraschi & Tavecchio 2009). 

Since the broad emission lines are produced by clouds photo- 
ionized by the radiation produced by the accretion disc, there is a 
direct relation between -Lblr and the accretion disc luminosity La. 
Therefore, measuring the broad line luminosities, we have infor- 
mation on the disc luminosity even when it is not directly visible, 
as often occurs in blazars whose optical continuum is dominated 
by the jet flux. In turn, by knowing Lblr and the bolometric jet 
luminosity, we can then study the relation between the jet and the 
accretion power. This in fact is the final aim of these studies. Earlier 
attempts to find the ratio between the jet and the accretion power 
were done by e.g. Celotti, Padovani & Ghisellini (1997) and by 
D'Elia, Padovani & Landt (2003): the novelty here is on one hand 
the way to estimate the jet power, and on the other hand the large 
number of sources for which the 7-ray detection ensures a good 
estimate of the jet power (or at least a good proxy for it), coupled 
with the large number of blazars present in the Sloan Digital Sky 
Survey (SDSS; York et al. 2000) with spectroscopic data. 

Another important ingredient for this line of research is the 
black hole mass, allowing to measure luminosities and powers in 
Eddington units. Besides allowing to compare objects with dif- 
ferent black hole masses, it allows to investigate if the accretion 
regime has indeed a transition, from radiatively efficient to ineffi- 
cient, when the mass accretion rate in Eddington units Al /M-^dd 
goes below a critical value (see e.g. Narayan, Garcia & McClin- 
tock 1997), and to see how this influences the jet power. For in- 
stance, the division between BL Lacs and FSRQs in the sam- 
ple of blazars detected during the first three months of the Fermi 
all sky survey (LBAS) seems to corresponds to disc luminosities 
Ld/Lsdd ~ 10~ 2 (Ghisellini, Maraschi & Tavecchio 2009; see 
also the earlier proposal about the division of FR I and FR II radio- 
galaxies in Ghisellini & Celotti 2001). 

For these reasons we are motivated to enlarge the original sam- 
ple of Gl 1, studying all blazars for which we can have information 
about their emission lines (as a proxy for the disc luminosity), their 
7-ray luminosity (as a proxy for the jet power), and their black hole 
mass. The two largest samples useful for this study are the SDSS 
and the Fermi 1LAC sample. In §2 we present the samples used 
for this work and in §3 we discuss how we have derived the broad 
line luminosities, or their upper limits. In §4 we present the rela- 
tion between Lblr and the 7-ray luminosities, and we discuss our 
findings in §5. 



2 THE SAMPLES 

We are interested in grouping a large number of blazars with reli- 
able measures of Broad Line Region and 7-ray luminosities. The 
Sloan Digital Sky Survey (SDSS), that provides the largest publicly 
available catalog of spectral objects, and the Large Area Telescope 
(LAT) onboard the Fermi Gamma-Ray Space Telescope constitute 
optimal devices to investigate in this direction. 

We tried to select for our analysis the largest group of blazars 
with reliable redshift and black hole mass measures. At first, we 
grouped a sample of optically selected quasars from the SDSS (sev- 
enth data release, DR7), that have been deeply analysed by Shen et 
al. (2011; hereafter Sll) and are Fenm-detected. Furthermore, in 
order to enlarge our analysis towards lower luminosities, we in- 
cluded in our study an optically selected group of BL Lac candi- 
dates. The BL Lacs that we took in consideration are supposed to 
be lineless, but their redshift and black hole masses have been de- 



rived in the work by Plotkin et al. (2011; hereafter PI 1) from the 
galaxy spectral absorption features. In the end, we tried to look for 
possible intermediate objects, that have been excluded both from 
SDSS DR7 quasar catalog and PI 1 BL Lacs catalog. For this pur- 
pose, we selected a small group of AGNs from the previous SDSS 
data release (DR6, Adelman-McCarthy et al. 2008) that are Fermi- 
detected. 

2.1 1LAC sample 

This is the AGN sample resulting after 1 1 months of operation of 
the Large Area Telescope (LAT) onboard the Fermi satellite (First 
LAT AGN Catalog, 1LAC, Abdo et al. 2010a). As revised in Ack- 
ermann et al. (2011), this sample is made by 671 sources at high 
Galactic latitude (|6j > 10°). The statistical significance of the 
Fermi detection was required to be TS>25 (TS stands for Test 
Statistics, see Mattox et al. 1996 for the definition. TS=25 approx- 
imately corresponds to 5<r). 

By requiring that the associations of the detected sources have 
a probability P ^ 80% and that there is only one AGN in the po- 
sitional error box of the Fermi detection, the Fermi team has con- 
structed the clean 1LAC sample, composed of 599 AGNs. Of these, 
248 are classified by the FermifLXT team as FSRQs and 275 as BL 
Lacs. The remaining sources are either of unknown blazar type (50) 
or non-blazar AGNs (26). We then focus on the 248 FSRQs and 
275 BL Lacs. All FSRQs have known redshift, while for about half 
of the BL Lac objects the redshift is still unknown. 

2.2 SDSS DR7 Quasar sample 

We have collected at first a group of optically selected quasars. The 
SDSS provides a quasar catalog out of the Data Release 7, that in- 
cludes 105783 objects and has been spectrally analysed by Sll. 
This sample includes quasars that have luminosities larger than 
Mi = -22 [i.e. I/L45100 A) = 10 44 erg s _1 ], that have at least 
one emission line with FWHM> 1000 km s _1 and a reliable spec- 
troscopical redshift. For this group of sources Sll calculated con- 
tinuum and emission line measurements around the Ha, H/3, Mgll 
and CIV regions, and derived virial black hole masses with differ- 
ent calibrations, along with what they consider the best estimate. 
The broad line luminosities and the best estimate of the black hole 
masses are considered for our work. 

From the cross-correlation with the 1LAC sample, we ob- 
tained a group of 49 Fermi detected and optically selected quasars. 
We excluded 3 objects because SI 1 do not provide a reliable black 
hole mass estimate. Therefore, the quasar sample under study in- 
cludes 46 objects, that are listed in TableQ] Note that requiring that 
the objects have a broad emission line with measurable FWHM au- 
tomatically excludes BL Lac objects with very weak or no emission 
lines. 

2.3 Plotkin et al. (2011) BL Lac sample 

We included in our analysis all the optically selected BL Lacs that 
are present in the work by PI 1. In their work, Pll start from an 
original sample of 723 BL Lac candidates with EW<5 A, that 
are included in the DR7 general catalog (for selection details, see 
Plotkin et al. 2010). From this original sample, they selected 143 
BL Lac candidates with reliable redshift limited to z < 0.4, that 
match to a FIRST and/or NVSS radio source and are radio-loud 
(FIRST stands for Faint Images of the Radio Sky at 20 cm, White 
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et al. 1997; NVSS stands for NRAO VLA Sky Survey, Condon et 
al. 1998). Then they apply a spectral decomposition in order to sep- 
arate the galaxy and AGN spectral components. Only 71 out of the 
143 BL Lacs could be successfully decomposed. For this smaller 
sample of 71 BL Lacs, black hole masses are derived from the M- 
cr* relation. 

A similar study on black hole masses of SDSS BL Lac objects 
has been performed in a work by Leon-Tavares et al. (2011). The 
authors started from a sample of BL Lacs included in the SDSS 
Data Release 5 (DR5) and radio-detected by FIRST. This original 
BL Lac sample was selected by Plotkin et al. (2008). Leon-Tavares 
et al. performed a spectral decomposition similar to the one in PI 1 
on objects in the redshift range 0.06 < z < 0.5. They obtain the 
black hole mass estimates for 78 BL Lacs, using the M-a, relation. 
The results of Leon-Tavares et al. (201 1) and those of PI 1 are very 
similar and the two works are consistent. We then decided to use 
the PI 1 data. 

We cross-correlated the Pll sample with the clean 1LAC 
sample. 10 BL Lacs out of the 71 are detected by Fermi, hence we 
have measures of their 7-ray luminosities. For the other 61 sources, 
we derived an upper limit on their 7-ray fluxes, based on the sen- 
sitivity limit of LAT for objects with T 7 ~ 2. Therefore, the upper 
limit in flux for these 61 BL Lacs is fixed at F p h = 5 x 10 -9 
ph cm~ 2 s _1 . The overall BL Lac sample under study includes 71 
objects, that are listed in Tablef2] 



2.4 SDSS DR6 sample 

We selected a last group of sources from the Data Release 6 (DR6) 
of the SDSS, in order to include in our analysis possible opti- 
cally intermediate objects, that may have been excluded from the 
Sll and Pll catalogs. Therefore, we cross-correlated the SDSS 
DR6 and the clean 1LAC and obtained a group of 20 additional 
sources, not contained in the samples mentioned above. Three of 
these sources have no reliable redshifts, and we excluded them 
from our sample. Since for this group of sources we do not have 
a black hole mass estimate, we have chosen to assign them an av- 
erage value of M — 5 x 1O 8 M0. This last group of intermediate 
blazars includes 17 objects, that are listed in Table[3] 

To sum up, the sample on which we work is composed by: 

• 46 Fermi-detected, optically selected FSRQs from Sll. These 
objects have detection both on Lblr and F 7 , and have black hole 
masses estimated by S 1 1 . 

• 10 Fermi-detected, optically selected BL Lacs from PI 1. Be- 
cause of their original selection, these sources do not show any 
emission line. Therefore, we calculated the upper limit on Fblr, 
while detections are available for L 7 . The 1LAC catalog provides 
for one of these sources an upper limit instead of a detection. PI 1 
provide mass estimates for all these objects. 

• 61 optically selected BL Lacs from Pll that are not Fermi- 
detected. As the other 10 sources from Pll, they have mass esti- 
mates (from Pll) and upper limits on Lblr, but in addition we 
calculated upper limits on L 7 , too. 

• 14 Fermi-detected objects that are included in the DR6 gen- 
eral sample. These objects do not show broad emission lines in their 
spectra, hence we calculated upper limits on their Fblr- In DR6 
there are no estimates of the black hole mass, hence we assigned to 
these objects an average mass value (M = 5 x 1O 8 M0). 

• 3 Fermi-detected blazars that are included in the DR6 sample. 
These objects show at least one emission line in their spectra, hence 



they have detections on both Fblr and F 7 . For one of them, we 
calculated a mass estimate from the FWHM of the H/3 line, while 
we assigned to the others an average mass value. 

In total, we have 49 objects with detections on both Lblr and L 7 , 
23 with upper limits on Lblr and detections on L 7 and 62 with up- 
per limits on both luminosities. In the following, when discussing 
the relation between the BLR and the 7-ray luminosity, we will 
add to our sample other 30 blazars studied in Gl 1, that are listed in 
TableU Of these 30 objects, 14 are FSRQs and 16 are BL Lacs (12 
LBLs and 4 HBLs). 29 have measured Fblr and 7-ray detections, 
while one has an upper limit on Fblr and a 7-ray detection. The 
total number of blazars with measured Fblr, F 7 and black hole 
mass is therefore 78. 



3 BROAD LINE LUMINOSITIES 

We have taken the luminosity of the emission lines of the blazars in 
the SDSS DR7 Quasar sample directly from the listed values in the 
Sll catalog. For calculating the total luminosity of the broad lines, 
we have followed Celotti, Padovani & Ghisellini (1997). Specifi- 
cally we set the Lya flux contribution to 100, the relative weight of 
Ha, U/3, Mgll and CIV lines respectively to 77, 22, 34 and 63 (see 
Francis et al. 1991). The total broad line flux is fixed at 555.76. The 
Fblr value or upper limit of each source has been derived using 
these proportions. In Tab.[T]we list these blazars reporting the type 
of lines used for calculating Fblr and the values of the estimated 
Fblr and the observed L 7 . When more than one line is present, 
we calculate the simple average of the Lblr estimated from each 
line. 

3.1 Upper limits on the broad line luminosity 

While the SDSS DR7 Quasar sample is selected in order to con- 
tain spectra with prominent broad emission lines, that have been 
measured by Sll, the other two samples include mostly lineless 
objects. In these cases, we need to derive upper limits on the line 
fluxes (ULF linc ). 

To this aim, the observed spectrum has been fitted with a 
power-law model, with the addition of possible absorption or nar- 
row emission features modeled as Gaussian profiles. Absorption 
lines are well visible in the spectra included in PI 1 BL Lac sample 
and the Pll work itself provides the variance of each absorption 
feature (a,). Such a modeling includes the lineless power-law and 
the narrow features. The broad emission line for which we want 
to obtain the upper limit is accounted as an additional Gaussian 
profile, with a variable flux value Fii nc and a FWHM fixed at the 
average value ufwhm = 4000 km s _1 , as suggested in Decarli et 
al. (201 1). This value is an average for all blazars, and it is consis- 
tent with the median FWHM values that can be obtained from the 
whole SDSS Quasar sample. Even though in the case of BL Lacs 
the average value is possibly slightly smaller, we prefer to maintain 
a larger average FWHM value, in order to derive more conservative 
(i.e. less stringent) upper limits. To define the ULf 11m we perform 
a x 2 test, varying the Fjj ne value until our model returns an unac- 
ceptable fit. Then we define the upper limit on the line flux as the 
Fii ne for which we obtain \ 2 > x 2 (99%). Over this critical value, 
the model is no more acceptable to fit the data, and we should ac- 
tually see a broad line emerging over the continuum, if present. In 
order to derive meaningful upper limits, we required a signal-to- 
noise ratio S/N > 5 in the wavelenght interval in which we per- 
formed our analysis. Hence we checked the signal-to-noise ratios 
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Table 1. Sources from the DR7 Quasar Catalog that are present in the 1LAC Fermi sample. Col. [1]: name; Col. [2]: right ascension; Col. [3]: declination; 
Col. [4]: redshift; Col. [5]: Logarithm of the black hole mass (in solar masses, best estimate from Sll); Col. [6]: lines measured by SI 1, from which -Lblr 
has been derived; Col. [7]: Broad Line Region luminosity (10 42 erg s _1 ), obtained from the line luminosities calculated by Sll; Col. [8]: 7-ray luminosity 
from Fermi data (10 45 erg s — 1 ), averaged on the first 1 1 months of Fermi operations. 



of the spectra in our sample, and we excluded B3 1432+422 (SDSS 
J143405.69+420316.0, from the DR6 sample), because its signal- 
to-noise ratio was S/N < 5 over the whole spectrum. Therefore, 
we are left with 16 DR6 objects, listed in Table[3] 

In principle, the process used to derive the upper limits could 
be applied to the four lines measured in the work by S 1 1 (i.e. Ha, 
H/3, Mgll and CIV). The objects included in the Pll sample all 
have z < 0.4, hence the ULFi ine can be derived only for the Ha 



and H/3 lines. We then derive the upper limits for these two lines. 
We applied the procedure also to the low redshift objects included 
in the DR6 sample. In 3 objects from this sample the redshift is 
sufficiently large to derive the upper limit on the Mgll line. 
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Name (SDSS J. . . ) 


z 


log M/Mq 


Lines 


UL Lblr 


ULL 7 












[10 42 ergs" 1 ] 


[10 44 ergs" 1 ] 


[10 44 ergs- 1 ] 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 



002200.95 +000657.9 





.306 


8.49 


Ho 


H 


3 


.65 


1.07 




005620.07 -093629.7 





.103 


9.01 


Ho 


II 


1 


.02 


0.95 




075437.07 +391047.7 





.096 


8.24 


Ho 


H/3 





.52 


0.81 




080018.79 +164557.1 





.309 


8.58 


Ho 


H/3 


5 


.68 


1.09 




082323.24 +152447.9 





.167 


8.80 


Ho 


H/3 


1 


.62 


2.69 




082814.20+415351.9 





.226 


8.83 


Ho 


II 


2 


.55 


5.24 




083417.58 +182501.6 





.336 


9.34 


H/3 




5 


.79 


13.18 




083548.14+151717.0 





.168 


7.94 


Ho 


H/3 


2 


.87 


2.75 




083918.74 +361856.1 





.335 


8.50 


H/3 




7 


.89 


13.18 




084712.93 +113350.2 





.198 


8.52 


Ho 


h/3 


3 


.38 




3.34 


085036.20 +345522.6 





.145 


8.61 


Ho 


Hri 


2 


.07 




1.74 


085729.78 +062725.0 





.338 


8.23 


H/3 




7 


.23 


13.48 




085749.80 +013530.3 





.281 


8.69 


Ho 


H/3 


5 


.85 


8.70 




090207.95 +454433.0 
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Ho 


H/3 


5 


.24 


9.33 




090314.70+405559.8 





.188 


8.28 


Ho 


Rd 


2 


.02 


3.54 




090953.28 +310603.1 





.272 


8.95 


Ho 


H/3 


6 


.54 


8.12 




091045.30 +254812.8 





.384 


8.51 


H/3 




9 


.92 


18.19 




091651.94+523828.3 





.190 


8.53 


Ho 


H/3 


3 


.07 


3.63 




093037.57 +495025.6 





.187 


8.48 


Ho 


H/3 


3 


.26 


3.46 




094022.44+614826.1 





.211 


8.57 


Ho 


H/3 


3 


.55 




11.57 


094542.23 +575747.7 





.229 


8.63 


Ho 


H/3 


3 


.23 




15.97 


101244.30 +422957.0 





.365 


8.67 


H/3 




13.30 


15.84 




102453.63 +233234.0 





.165 


7.46 


Ho 


Hri 


2 


.68 


2.63 




102523.04 +040228.9 





.208 


8.18 


Ho 


H/3 


2 


.74 


4.36 




103317.94+422236.3 
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8.59 


Ho 


II 


3 


.18 


4.57 
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.304 


8.70 


Ho 


H/3 


6 


.45 


10.47 




104149.15 +390119.5 





.208 


8.55 


Ho 


H/3 


3 


.23 


4.36 




104255.44+151314.9 





.307 


7.81 


Ho 


H/3 


4 


.56 


10.71 




105344.12+492955.9 





.140 


8.47 


Ho 


H/3 


1 


.68 




3.03 


105538.62 +305251.0 





.243 


8.43 


Ho 


H/3 


4 


.10 


6.30 




105606.61 +025213.4 





.236 


8.11 


Ho 


H/3 


3 


.43 


5.88 




105723.09 +230318.7 





.378 


8.32 


H/3 




9 


.91 


17.37 




112059.74 +014456.9 





.368 


9.60 


H/3 




9 


.81 


16.21 




113630.09+673704.3 





.134 


8.30 


Ho 


H/3 


1 


.07 




2.03 


114023.48+152809.7 





.244 


9.46 


Ho 


H/3 


4 


.43 


6.30 




114535.10-034001.4 





.168 


8.27 


Ho 


H/3 


2 


.08 


2.75 




115404.55-001009.8 





.254 


8.36 


Ho 


H/3 


3 


.54 




5.53 


115709.53 +282200.7 





.300 


9.20 


Ho 


H/3 


5 


.60 


10.23 





Table 2. BL Lacs from the work by Pll. Col. [1]: SDSS name; Col. [2]: redshift; Col. [3]: Logarithm of the black hole mass (Pll, in solar masses); Col. [4]: 
lines from which we derived the upper limits, as described in Section [3~T1 Col. [5]: upper limit on the Broad Line Region luminosity (10 42 erg s -1 ), obtained 
from the UL on line fluxes; Col. [6]: upper limit on the 7-ray luminosity obtained from the Fermf'-LAT sensitivity limit, in units of 10 44 erg s -1 ; Col. [7]: 
7-ray luminosity from Fermi data (10 44 erg s — 1 ), averaged on the first 1 1 months of Fermi operations. 



3.2 The distribution of black hole masses 

As a byproduct of our study, we have collected (from the S 1 1 and 
Pll samples) a large number of estimates for the black hole masses 
both in FSRQs and BL Lacs. FSRQs show a distribution skewed to- 
wards larger masses than BL Lacs: (log A/fsrq) = 8.88 ± 0.40; 
(log Mbl Lac) = 8.57 ± 0.37, while the average of all masses is 
(log M a n) = 8.70 ±0.41. We believe that at least in part this is due 
to a selection effect, since most BL Lacs come from the Pll sam- 
ple, therefore they have been selected to be at z < 0.4. Assuming 
that very large black hole masses are rarer than smaller ones, one 
expects that the largest masses are found only when considering 
relatively large redshifts. This is illustrated by Fig. [T] that shows 
the black hole masses as a function of redshift for BL Lacs and FS- 
RQs. Apart from few exceptions, the BL Lacs extend out to z ~ 0.4 



(by construction, given the redshift limit of the Pll sample), while 
FSRQs cluster around z ~ 1. 



4 THE L B lr-£ 7 RELATION 

Fig.|2]is the key result of our work. It shows the luminosity of the 
BLR as a function of the observed 7-ray luminosity, both mea- 
sured in Eddington units. Arrows correspond to upper limits. Dif- 
ferent symbols correspond to blazars belonging to different sam- 
ples, as labelled. Note that we have also added the blazars studied 
in Gil, but omitting the objects in common. Fig. [2] shows a clear 
trend. Since the range in black hole masses is relatively narrow, we 
obtain a similar trend when plotting Lblr vs L t . We have quan- 
tified it first by using the Kendall non-parametric test considering 
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Table 2. BL Lacs from the work by PI 1. *: the upper limit is from the Abdo et al. (2010a) list. Col. [1]: SDSS name; Col. [2]: redshift; Col. [3]: Logarithm 
of the black hole mass (PI 1, in solar masses); Col. [4]: lines from which we derived the upper limits, as described in Section |3~71 Col. [5]: upper limit on 
the Broad Line Region luminosity (10 42 erg s -1 ), obtained from the UL on line fluxes; Col. [6]: upper limit on the 7-ray luminosity obtained from the 
Fermi-LAT sensitivity limit, in units of 10 44 erg s — 1 ; Col. [7]: 7-ray luminosity from Fermi data (10 44 erg s _1 ), averaged on the first 1 1 months of Fermi 
operations. 



the detected sources (i.e. excluding upper limits). The Kendall r in 
this case is listed in Tab. [5] The correlation is significant both con- 
sidering log I/blr vs log L 7 and when measuring these quantities 
in Eddington units. Then we have considered the common depen- 
dence upon the redshifts of both luminosities, and applied the par- 
tial Kendall correlation analysis as described in Akritas & Siebert 
(1996). The correlation is still significant, although with a smaller 
r. We then included the upper limits, and repeating the same anal- 
ysis we verify that the value of r is now greater. 

Finally, we applied a simple least square fit and performed a 
partial correlation analysis (see Eq. 1 of Padovani 1992), account- 
ing for the common dependence on the redshift and on the black 
hole mass of the plotted quantities. The results are listed in Tab. [5] 
In this case we have excluded all upper limits from the analysis. 

Before discussing the implications of this correlation, there are 
a few caveats to remind, concerning possible important selection 
effects: 

• In the 1LAC catalog there are many detected sources without 



a known redshift. As discussed in Abdo et al. (2010b) and in Gl 1, 
if these sources will turn out to be at 2 ~ 2, then their 7-ray lu- 
minosities would be huge. If the absence of broad emission lines is 
due to their intrinsic weakness, then these blazars would be located 
in the bottom right part of Fig. [2] and they would be clear outliers 
of the found correlation. If, instead, the absence of lines is due to 
a particularly strong non thermal continuum, then Lblr could be 
large, locating these objects in the "FSRQs quadrant". 

• We have clear examples of blazars varying their 7-ray lumi- 
nosity by more than 2 orders of magnitude. It is very likely that the 
present samples of 7-ray detected blazars preferentially include ob- 
jects in their high statqj. The 7-ray luminosity we have considered 
is the average over 1 1 months (therefore the short term variability is 



1 This would also explain why the radio and 7-ray fluxes are correlated, 
even if only a relatively small fraction of radio-loud AGN with a flat spec- 
trum are detected in 7-rays; see e.g. Ghirlanda et al. (201 1). 
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Table 3. Sources from the DR6 Catalog that are present in the 1LAC Fermi sample. Col. [1]: name; Col. [2]: right ascension; Col. [3]: declination; Col. [4]: 
redshift; Col. [5]: lines measured or for which we derived the upper limits as described in Section [3~71 Col. [6]: upper limit on the Broad Line Region luminosity 
(10 42 erg s — 1 ), obtained from the UL on line fluxes; Col. [7]: luminosity of the Broad Line region (10 42 erg s _1 ); Col. [8]: 7-ray luminosity from Fermi data 
(10 45 erg s- 1 ), averaged on the first 1 1 months of Fermi operations. The black masses are not available, hence a medium mass value (M = 5 X 10 8 M o ) 
has been assigned to all of them. a : PKS 1 106+023 has the H/3 line measured in DR6 Catalog, hence its mass can be estimated by the Chiaberge & Marconi 
relation (2011): a value of M = 4 X 1O 7 M is obtained. 
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Figure 1. Black hole masses as a function of redshift. Empty circles are 
BL Lacs, filled circles are FSRQs. BL Lacs of our sample have significantly 
smaller redshift than FSRQs. 

averaged for), but blazars can indeed be variable over longer peri- 
ods. This variability introduces an inevitable dispersion around the 
correlation line. 

• Misaligned jets should be weaker 7-ray sources than their 
aligned counterpart, but they would show the same emission line 
luminosities. Therefore weak 7-ray sources must exist, populating 
the region to the left of the interpolating line of Fig. [2] However, 
these sources would be classified as radio-galaxies (see Abdo et al. 
2010c), and not aligned blazars, although some overlap might exist. 

Despite the presence of the above caveats, the apparent cor- 



relation between the BLR and the jet 7-ray luminosity is certainly 
intriguing, since it would prove the importance of the emission line 
photons in the production of high energy 7-rays and, more impor- 
tantly, it would point towards a relation between the accretion rate 
and the jet power. This relation is not direct, however, since the ob- 
served 7-ray luminosity can be considered a rather poor proxy of 
the jet power and the disc luminosity is linearly related to the ac- 
cretion rate only for "standard" optically thick geometrically thin 
accretion disc. This will be discussed more in the next section. 

Fig. [2] shows that BL Lac objects are rather neatly divided 
from FSRQs: with few exceptions, all FSRQs have Lblr/ ^Edd > 
5 x 10~ 4 and all BL Lacs are below this value. The correspond- 
ing dividing 7-ray luminosity is L 7 /LEdd ~ 0.1. We derived this 
apparent "divide" considering only the sources for which we have 
a detection of the BLR and the 7-ray luminosities. In other words, 
we first excluded the upper limits from the analysis. The paucity 
of data does not allow us to reach a firm conclusion about the ex- 
act value of this divide, but, reassuringly, when we include all the 
upper limits they lie in the "correct" quadrant of the plane. 

On the other hand, most upper limits correspond to BL Lacs 
with z < 0.4, so an issue remains: the divide between BL Lacs 
and FSRQs could be partly due to a segregation in redshift, if all 
BL Lacs are at low redshift and FSRQs at high redshift. To verify 
this, we plot in Fig.[3]the BLR luminosity (in Eddington units) as a 
function of redshift. It can be seen that there are detected BL Lacs 
with LBLR/^Edd < 5 x 10~ 4 at relatively large redshifts. This 
is a hint that the divide is real, but, again, the paucity of points 
precludes a firmer conclusion about this possible selection effect. 
Moreover, most of the upper limits come from the PI 1 sample, that 
by construction selects only BL Lacs at z < 0.4. This limit in 
redshift possibly introduces a bias in the dividing value, as can be 
seen in Fig|3] Nevertheless, we reiterate that the upper limits did not 
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Table 4. Blazars from Gil. Col. [1]: name; Col. [2]: right ascension; Col. [3]: declination; Col. [4]: redshift; Col. [5]: Logarithm of the black hole mass (in 
solar masses); Col. [6]: Broad Line Region luminosity (10 42 erg s — 1 ); Col. [7]: 7-ray luminosity (10 45 erg s — 1 ). 



take part in the determination of the divide, hence this bias does not 
completely compromise the result. 

If real, the found divide would be in agreement with what 
found studying the distribution of bright Fermi detected blazars in 
the 7-ray spectral index - 7-ray luminosity plane (for them Ghis- 
ellini Maraschi & Tavecchio 2009 proposed a "divide" between BL 
Lacs and FSRQs around I/d/LEdd ~ 10 _ ), and to what more re- 
cently found by Gil using a sample of bright blazars much more 
limited in number than what we use here. The value of the divide 
found here would also be consistent with the division between FR 
I and FR II radio-galaxies, found with a completely different ap- 
proach by Ghisellini & Celotti (2001). 

Fig. [2] shows also that all blazars form a continuous family, 
with no apparent "discontinuity" (or sign of bimodality). Exclud- 
ing upper limits, we find that Lblr oc L 7 (normalizing or not to 
Eddington). One can then wonder if it is still meaningful to divide 
BL Lac objects from FSRQs, since in Fig.|2]they form a continuous 
distribution. In other words: are BL Lacs and FSRQs characterized 



by some different fundamental properties, or are they simply at two 
sides of a continuous distribution of properties? An example can il- 
lustrate this point: suppose, as suggested by Ghisellini, Maraschi & 
Tavecchio (2009) that the accretion disc in BL Lacs is radiatively 
inefficient, while it is efficient in FSRQs. This is a fundamental dif- 
ferent property, although it concerns the accretion disc, not the jet. 
Another example: suppose that jets in BL Lacs are made by pure 
electron-positron plasmas, while the jets in FSRQs are made by 
normal electron-proton plasmas. This, too, should be considered 
a fundamental different property. If instead all blazars have radia- 
tively efficient accretion disc and their jets are all made by electrons 
and protons, then they look different only because they have differ- 
ent overall powers, and this in turn might also explain why their 
SED is different, without the need of anything fundamental divid- 
ing them. 

The discussion below is dedicated to this issue, focussing in 
particular to the proposed "divide" of BL Lacs and FSRQs in terms 
of the mass accretion rate in Eddington units. 
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Figure 2. Luminosity of the broad line region (in Eddington units) for the sources from our samples and from Gil as a function of the 7-ray luminosity (in 
Eddington units). Different symbols correspond to different samples or a different classification of the sources, as labelled. The three (violet) asterisks are the 
only sources with visible broad emission lines from the DR6 Fermi detected sample. The three labelled triangles have their synchrotron emission dominating 
over the thermal emission in their SEDs. Hence, to avoid black hole mass estimate errors, possibly occurred in the SI 1 automatic calculation, we also assigned 
them an average Mbh value (Mbh = 5 x 10 s Mq). These changes are highlighted by the thick (black) segments ending to the black triangles (corresponding 
to the average A/bh value). The grey stripe indicates the luminosity "divide" between FSRQs and BL Lacs at -^BLR/^Edd ~ 5 x 10 — 4 . 



5 DISCUSSION 

Our study concerns Lblr, L 7 , and the black hole mass. In the fol- 
lowing we will discuss how we can use Lblr to find the disc lumi- 
nosity L<j, and how to use L 7 to find a proxy for the jet power Pj ot 
and the mass accretion rate M. The black hole mass is of course 
used to normalize all powers to the Eddington luminosity. 

Lblr — >■ Ld — > M — For radiatively efficient accretion 



discs, the BLR luminosity is a direct measure of the disc lumi- 
nosity Ld, since, on average, Ld ~ 10 Lblr (see e.g. Baldwin & 
Netzer 1978; Smith et al. 1981). Radiatively efficient (i.e. Shakura- 
Sunjaev 1973) disc should occur for rh = Af/A/Edd > rn c . Defin- 
ing A/sdd = LEdd/c 2 (without the efficiency factor), then m c 
should be close to 0.1 (Narayan & Yi, 1995). Another hypothe- 
sis suggests the lower value m c ~ 10 -4 for the radiative transition 
(Sharma et al. 2007). If the disc emits as a black-body at all radii, 
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Table 5. The first part of the table reports the results of the non-parametric Kendall test for the complete sample, taking into consideration at first only the 
sources with both Lblr an d L 7 detected, then including upper limits. We list also the results when accounting for the common dependence on redshift. The 
bottom part of the table reports the results of a partial correlation analysis using a least square fit. We have excluded upper limits from the analysis. Correlations 
are of the form log y = m log x + q. The listed slopes m refer to the bisector (of the two correlations x vs y and y vs x). N is the total number of objects, r 
is the correlation coefficient obtained from the analysis. P is the probability that the correlation is random. 
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Figure 3. The broad line luminosity (in Eddington units) as a function of 
redshift. Same symbols as in Fig.|2] 



then most of the power is emitted in the far UV, and we can approx- 
imate the photo-ionizing luminosity with the entire L d . 

When rh < rh c , the disc should become radiatively ineffi- 
cient, because the particle density of the accretion flow becomes 
small, and the energy exchange timescale between protons and 
electrons becomes smaller than the accretion time. If this occurs, 
the disc bolometric luminosity decreases. Narayan, Garcia & Mc- 
Clintock (1997) proposed that in this regime Ld oc M 2 . In this case 
the disc becomes hot, inflates, and it does not emit black-body ra- 
diation. As a consequence, Li on <H Ld- According to Mahadevan 
(1997, see his Fig. 1), the decreasing fraction of the ionizing lu- 
minosity is as important as the decrease of the overall efficiency 77 
(defined as Ld = r/Mc 2 ). In the example shown by Fig. 1 of Ma- 
hadevan (1997), Li on oc M 3 . If this were true, we expect also the 
broad emission line luminosity to have the same dependence on M 
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Figure 4. The 7-ray luminosity in the LAT band as a function of the jet 
power P r (both in Eddington units). The latter is the power that the jet 
has spent to produce the (bolometric) radiation we see, and is given by 
P r ~ ibol/F 2 . To derive it we have used blazars detected by Fermi for 
which we have constructed the SED and estimated the bulk Lorentz factor 
(G10; Gil; Tavecchio et al. 2010). The grey solid line shows the result 
(bisector) of a least square fit. 



when Al goes sub-critical. Note, however, that the SED calculated 
by Mahadevan (1997) could be largely affected by the presence 
of extras sources of seed photons for the thermal Componization, 
besides the assumed cyclotron-synchrotron emission. These extra 
seed photons (for instance coming for some cool part of the disc) 
could enhance the UV emission both by enhancing the Compton 
scattering and by cooling the hot emitting electrons. So we regard 
the Lblr oc Lion oc M ' relation as an indication, but without 
excluding other possibilities. In practice, for m < m c , we will 
consider both Lblr oc M 2 and Lblr oc M ' . 
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L 7 — > Pj c t — > M — A lower limit on the jet power Pj e t is 
P ^ P ■> Lh ° l 

Jjct ^ r-c ^ 2 

where F is the bulk Lorentz factor (see Celotti & Ghisellini 2008 
and Ghisellini & Tavecchio 2009), Lboi is the jet bolometric lumi- 
nosity, measured assuming isotropic emission. P r is the power that 
the jet has to spend to produce the radiation we see. If Pj ct ~ P r , 
then the jet uses its entire power to produce radiation, including 
its kinetic power, and it should stop. Fitting the SED of a large 
number of blazars detected by Fermi with a simple one-zone lep- 
tonic model (as the one used in the Appendix) returns values of 
r contained in a narrow range, around 13-15 (see Ghisellini et al. 
2010; hereafter G10), consistent with the values one derives from 
the superluminal motion. The same model also yields the number of 
emitting electrons required to account for the radiation we see, and 
by assuming one proton per electron, it yields Pj et ~30-100P r , on 
average. We can conclude that Pj c t is robustly bound to be larger 
than P r ~ Lboi/200 and it can be a factor 30-100 larger than that. 

The key result obtained in G10 (confirming earlier results in 
Celotti & Ghisellini 2008 and then also confirmed by Gil), is that 
Pjct ~ Ale 2 . Since P r is linearly related to Pj e t, and it is a model- 
independent quantity, we can safely use P r as a proxy for Pj ot and 
M. To measure P r ~ iboi/T 2 we should construct the SED of 
all our blazars. On the other hand, we can take advantage from the 
sample already studied in the literature, to find if there is a robust 
relation between P r and the 7-ray luminosity L 7 in the Fermi/LAT 
energy range, that is the quantity most readily available to us. Such 
a robust relation indeed exists, and it is shown in Fig. [4] Fitting 
it with a simple least square method and taking the bisector we 
obtain: 



Log 
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Note that the two quantities are not linearly related. At low values 
of L 7 /I/Edd, the 7-ray luminosity underestimate P r (either nor- 
malizing or not to the Eddington luminosity). As a consequence, 
the L 7 /Z/Edd values for low luminosity BL Lacs underestimate the 
jet power, and in turn the mass accretion rate. This occurs because 
L 7 is not a good indicator of Lboi for low luminosity BL Lacs, that 
have their high energy SED peaking in the TeV band, and that have 
an important (often dominant) synchrotron component. 

i-BLR/^Edd vs Pr/I/Edd — Since the above correlation ap- 
pears rather robust, we can use it to obtain P r /Z/Edd from our 
values of I/ 7 and black hole mass, without the need to construct 
the SED and modelling it. In turn, we can think at the obtained 
P-/£Edd values as proportional to the mass accretion rate in Ed- 
dington units. This is done in Fig. [5] where we plot Lblr/L-e^ as 
a function of P-/I/Edd- 

For m > m c , we expect that the BLR luminosity, in Edding- 
ton units, depends linearly on the normalized accretion rate, and 
therefore L B Ln/L Edd OC 771 OC Pr/^Edd- 

Below 771 ~ 77i c , we expect that the disc becomes radia- 
tively inefficient. The disc bolometric luminosity becomes propor- 
tional to 771 2 , while the ionizing luminosity can follow an even 
steeper relation with 771 (i.e. oc m 3 ' 5 , see above). The grey stripe 
at small I/BLR/I<Edd is proportional to P r /-I/Edd in the top panel, 
to (Pr/^Edd) 2 in the mid panel and to (Pr/P/Edd) 3 ' 5 in the bot- 
tom panel. We mainly consider the Narayan & Yi hypothesis of 
m c ~ 0.1. The break value has been chosen to correspond to 
Lblb./ L-Ea ~ 5 x 10 -4 , and allowing for considerable scatter 
around this value. Note that, if we take into account the m c ~ 10 -4 
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Figure 5. Luminosity of the broad line region (in Eddington units) for the 
sources from our samples and from Gi l as a function of P r (in Eddington 
units). Symbols are the same as Fig. [2] but without upper limits and with 
grey stripes superimposed. In all three panels we assume that P r tracks M. 
If this is true, we can see how the broad line luminosity is related to rh = 
Af/Afsdd, or, equivalently, to Pr/^Edd- At large values of I/BLR/^Edd 
the grey stripe of all panels corresponds to ^blr/ ^Edd °£ Pr/^Edd- hi 
the top panel the grey stripe continues to show a linear relation also at low 
values of Lblr/ -^Edd, while in the mid panel the grey stripe becomes 
quadratic below some critical value (here L^/L^dd = 0-1 is assumed, 
namely the value dividing BL Lacs from FSRQs). In the bottom panel, the 
grey stripe becomes I/BLR/^Edd oc (-Pr/I/Edd) 3 ' 5 at l° w values, to ac- 
count for the deficit of ionizing UV photons in radiatively inefficient disc. 



hypothesis (Sharma et al. 2007), the break values would be outside 
the luminosity range of our sample. 

Fig. [5] shows that we cannot yet distinguish among the differ- 
ent cases, even if there is some preference for the Lblr/ L^dd oc 
(Pr/LEdd) 2 case. 

There are BL Lacs of the HBL type, such as Mkn 421, Mkn 
501 and 2005^189, that do show broad emission lines (and the pro- 
totypical BL Lac object too, i.e. BL Lac itself), with Lblr around 
10 42 erg s _1 and I/BLR/^Edd around 10 -5 . These are the objects 
at the low (bottom-left) end of Fig. [2] Fig. [4] and Fig. [5] We need 
more low power BL Lacs to investigate if they are indeed associ- 
ated with radiatively inefficient discs. 

Note that this is not required by the basic explanation of the 
blazar sequence and by the proposed luminosity "divide" between 
BL Lacs and FSRQs. In fact, the change of the observed SED along 
the blazar sequence, interpreted as a radiatively cooling sequence 
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(Ghisellini et al. 1998), requires that in BL Lacs the emission lines 
are not important as seed photons for the inverse Compton process. 
This can be the case even if the lines are present, if the dissipation 
region occurs outside Rblr: in this case the BLR photons are seen 
in the comoving region red-shifted and time-diluted, and the EC 
process can be negligible. The relation between the size of the BLR 
and the ionizing luminosity ensures that in BL Lacs the size of the 
BLR is much smaller than in FSRQs, even if the black hole mass is 
similar. If dissipation occurs always at 7? diaa ~ 10 3 Schwarzschild 
radii, then in BL Lacs we easily have Rdiss > Rblr. (and emit- 
ting lines are negligible for the formation of the high energy spec- 
trum), while in more powerful FSRQs we have -Rdiss < Rblr, 
with a corresponding enhancement of the EC process. Earlier works 
(Celotti & Ghisellini 2008; G10, Gil) have shown that i? d iss is of 
the order of ~ 10 3 Schwarzschild radii in all objects. Requiring 
that the size of the BLR is a factor / smaller than this, and using 
-Rblr = 10 17 L^/% 5 cm, we obtain 

Rblr < /Rdiss ^ < 6.9 x 10- 3 / 2 M 8 f^=-V 

(3) 

where M — 10 8 M 8 solar masses. We obtain, in this case, a value 
for the divide in agremeent with the observed I/BLR/LEdd ~ 
5 x 10~ 4 (for Lblr ~ O.ILd and / smaller than, but close to 
unity), but dependent on the black hole mass. The dependence on 
the black hole mass would produce some blur in the division, that 
is not inconsistent with what we see. 



6 CONCLUSIONS 

In this work we have studied those blazars detected by Fermi/LAT 
and present in the SDSS optical survey, for which the redshift is 
known and there is a black hole mass estimate. From the broad 
emission line luminosities (or their upper limits) we have calculated 
the luminosity of the entire Broad Line Region, used as a proxy for 
the luminosity of the accretions disc. We could find values for both 
the BLR and the 7-ray luminosity for 78 blazars, values for Lblr. 
and upper limits on L 7 for 23 blazars, and upper limits on both 
the quantities for 62 sources. Our results can be summarized as 
follows: 

(i) The luminosity of the BLR correlates well with the 7-ray lu- 
minosity in the Fermi/LAT energy range. The correlation is linear, 
irrespective if the above luminosities are normalized to Eddington 
or not. All upper limits (not used to find the correlation) are consis- 
tent with the correlation itself. 

(ii) BL Lac objects and FSRQs occupy different regions 
of the LBLR/L E dd-L 7 /L E dd plane, with a division at about 
Lb-lb, / L Edd ~5x 10~ 4 . This confirms, with an enlarged sample, 
earlier results. Nevertheless, since the sample is still rather poor of 
sources with detections on both Lblr and L 7 and instead rich of 
upper limits, this "divide" still needs further studies with a more 
populated sample. 

(iii) For objects (analyzed in previous works) of known L 7 , P r , 
and black hole mass, there is a strong correlation between the two 
quantities, both using absolute values and normalizing them to the 
Eddington luminosity: (L 7 /L E dd) oc (Pr/I/Edd) 1 ' 28 - As a conse- 
quence, the 7-ray luminosity (in the Fermi/LAT energy range) can 
be used to estimate P r which is a robust proxy for the jet power. 

(iv) The relation between the strength of the emission lines and 
the accretion rate can be used to test radiatively inefficient disc 
models and the prediction about the production, in these discs, of 



the ionizing luminosity. Our results are too primitive to draw strong 
conclusions, but there is the possibility that at low accretion rates 
the produced ionizing UV luminosity is larger than expected. 

(v) The division between BL Lacs and FSRQs could be due to 
the transition between a radiatively inefficient disc to a standard 
(Shakura-Sunyaev) disc. Alternatively, it can be due to the dissipa- 
tion region of the jet being located outside or inside the BLR. 
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7 APPENDIX 

7.1 Spectral energy distribution 

We have characterized the Spectral Energy Distribution (SED) of 
the six sources for which we have both the spectroscopic optical 
data and the detection by Fermi. 

To this aim we have collected the data from the NASA Ex- 
tragalactic Database (NED) and including the LB AS and 1LAC 
Fermi/LAJ data (Abdo et al. 2010a; 2010b). 

7.2 The model 

To model the SED we have used the leptonic, one-zone syn- 
chrotron and inverse Compton model, fully discussed in Ghisellini 
& Tavecchio (2009). 

In brief, we assume that in a spherical region of radius R, 
located at a distance i?di ss from the central black hole, relativistic 
electrons are injected at a rate Q(7) [cm -3 s _1 ] for a finite time 
equal to the light crossing time R/c. For the shape of Q(-y) we 
adopt a smoothly broken power law, with a break at 7b: 



0(7) 



(7/7b)" 



1 + (7/7b)~ sl+s 



(4) 



The emitting region is moving with a velocity j3c corresponding to 
a bulk Lorentz factor T. We observe the source at the viewing angle 
8 V and the Doppler factor is 5 — 1/[T(1 — /3cos# v )]. The mag- 
netic field B is tangled and uniform throughout the emitting region. 
We take into account several sources of radiation externally to the 
jet: i) the broad line photons, assumed to re-emit 10% of the ac- 
cretion luminosity from a shell-like distribution of clouds located 
at a distance Rblh = W 17 L\ % 5 cm; ii) the IR emission from 

a dusty torus, located at a distance Rm = 2.5 x 10 18 L^% 5 cm; 
iii) the direct emission from the accretion disc, including its X-ray 
corona; iv) the starlight contribution from the inner region of the 
host galaxy; v) the cosmic background radiation. All these contri- 
butions are evaluated in the blob comoving frame, where we cal- 
culate the corresponding inverse Compton radiation from all these 
contributions, and then transform into the observer frame. The lat- 
ter two contributions are negligible for our sources. 

We calculate the energy distribution AT (7) [cm -3 ] of the emit- 
ting particles at the particular time R/c, when the injection process 
ends. Our numerical code solves the continuity equation which in- 
cludes injection, radiative cooling and e pair production and re- 
processing. Our is not a time dependent code: we give a "snapshot" 
of the predicted SED at the time R/c, when the particle distribution 
N(y) and consequently the produced flux are at their maximum. 

To calculate the flux produced by the accretion disc, we adopt 
a standard Shakura & Sunyaev (1973) disc (see Ghisellini & Tavec- 
chio 2009). 

The resulting SEDs and models are shown in Fig.[6] Fig.UJand 
Fig. [8] and the model parameters are reported in Tab.|6]and Tab. [7] 



7.3 Specific sources 

We here discuss briefly a few objects for which we do have both the 
information on the BLR and the 7-ray luminosity. These sources 
lye close to the "intermediate zone" between BL Lacs and FSRQs. 
We have constructed their SED and modeled it through a simple 
one-zone leptonic model, as described in Section 17711 in order to 
classify them. To this aim we adopt the same SED-based clas- 
sification scheme discussed in Gil and originally introduced by 
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Padovani & Giommi (1995). In brief, we can classify the object as 
a FSRQ if the 7-ray luminosity is dominating the electromagnetic 
output and if the X-ray spectrum is flat (X-ray energy spectral in- 
dex dj < 1); it is a Low frequency peaked BL Lac (LBL) if the 7- 
ray luminosity is comparable to the synchrotron one and if ct x < 1, 
and high frequency peaked BL Lac (HBL) if the 7-ray luminosity 
is comparable or less than the synchrotron one and a x > 1. Fig. [6] 
Fig.[7]and Fig.[8]show the SED of these sources. We here summa- 
rize our finding. 

0937+5008 — This source is included in the Sll catalog and it 
is detected by Fermi. Looking at the SED it can be clearly classi- 
fied as a FSRQ, as suggested also from the evident broad Ha and 
H/3 lines visible in the SDSS spectrum. From the SED modeling, 
a synchrotron contamination of the optical continuum is visible. 
Therefore, the automatic virial black hole mass estimate performed 
by S 1 1 can be imprecise. Hence, we chose to assign to its black 
hole an average mass value M = 5 X 10 8 Mq. 

1040+24 — This is a DR6 Fermi detected source, hence the black 
hole mass is not measured. We can assume an average value of 
M = 5 x 1O 8 M . In the SDSS spectrum, a broad Mgll line is 
clearly visible. From the SED modeling, this source can be classi- 
fied as a LBL. This means that the thermal continuum can be highly 
contaminated by the synchrotron emissions, although some broad 
emission lines clearly emerge. Moreover, the disc seems to be only 
partially covered, and the synchrotron component looks really vari- 
able. This could result in a variable line EW. 

1055+01 — This source belongs to the DR6+1LAC sample, hence 
the mass is not measured and we assume the average value M = 
5 x 10 8 M Q . As in the case of 1040+24, a broad Mgll line is vis- 
ible, but it is narrower than the usual broad emission line width 
(FWHM~ 2500 km s" 1 ). This might suggest a small black hole 
mass, or maybe the line is partially covered by the continuum and 
hence the measure is uncertain. The SED, indeed, shows that the 
accretion disc contribution is mostly covered by the synchrotron 
emission. Overall, we can classify this source as a FSRQ, even if 
the disc emission is dominated by the synchrotron radiation. 

1106+023 — It is a DR6 source detected by Fermi. In this case 
we can estimate its mass with the Chiaberge & Marconi (2011) 
relation. In fact the FWHM of its broad H/3 line is reported in the 
DR6 catalog. The mass estimate that we obtain is very small (M = 
4 x lO 7 Af0). From the SED, we can see that the disc emerges from 
the synchrotron component, hence the emission lines can be seen 
clearly. From the SED we can classify this source as a LBL. 

1117+2014 — This source is present in the Sll catalog, hence it 
is supposed to be a quasar. However, in the SDSS spectrum, the 
H/3 line is really faint. Indeed, from the SED modeling, we can 
classify this source as an HBL. The SED shows clearly that the 
synchrotron component largely dominates the disc emission, hence 
the automatic mass estimate performed by SI 1 (M = 4 x 10 8 Mq) 
is not accurate. We then assume for this source a value of the black 
hole of equal to the average value. 

1509+022 — This is a source included in the Sll catalog. Look- 
ing at the SED, we can classify it as a FSRQ. However, also in 
this case the continuum appears contaminated by the synchrotron 
component. Moreover, the Sll results about this source are quite 
unclear. The equivalent widths reported in S 1 1 (EW~ 25 — 27 A), 
do not seem to be recognizable in the spectrum (the lines are hardly 
visible). Therefore, we think that the Sll black hole mass estimate 
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Figure 6. SED of 0937+5508 and 1040+24 and the fitting model. The dot- 
ted line corresponds to emission from the accretion disc, the IR torus (if 
present) and the X-ray corona, the thin (green) solid line is the synchrotron 
component, the long dashed line the SSC emission, the dot-dashed line the 
EC contribution. The thick (blue) solid line) is the sum. 
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Name 


z 


Pdiss 


M 


PbLR 


P.' 


L d 


B 


r 


e v 


7b 


7max 


si 


S2 


r 1 1 
L 1 J 


i z j 




L^J 






I ' 1 


t°j 


run 
[ >\ 


r i oi 


i 1 1 j 




n ^i 
i 1 -' i 


ri4i 


0937+5008 


0.275 


90 (600) 


5e8* 


34.6 


2e-3 


0.12(1.6e-3) 


0.14 


14 


3 


400 


5e3 





2.5 


1040+23 


0.56 


97.5 (650) 


5e8* 


51.2 


1.5e-3 


0.26 (3.5e-3) 


0.67 


13 


3 


100 


9e3 





2 


1055+01 


0.888 


105 (700) 


5e8* 


98.7 


0.02 


0.98 (0.013) 


3.1 


11 


3.7 


400 


9e3 


1 


2.5 


1106+023 


0.157 


43.2 (3.6e3) 


4e7 


21.9 


3e-3 


0.048 (8e-3) 


0.22 


13 


3 


6e3 


6e3 


2 


2 


1117+2014 


0.138 


105 (700) 


5e8* 


13.7 


5.e-5 


0.019 (2.5e^l) 


0.7 


15 


2 


3e4 


1.5e5 


0.5 


3 


1509+022 


0.2194 


120 (800) 


5e8* 


34.6 


0.022 


0.12(1.6e-3) 


0.13 


11 


6 


100 


2e4 





2.3 



Table 6. List of parameters used to construct the theoretical SED. Col. [1]: name; Col. [2]: redshift; Col. [3]: dissipation radius in units of 10 15 cm and (in 
parenthesis) in units of Schwarzschild radii; Col. [4]: black hole mass in solar masses, the asterisk means that the mass is assumed; Col. [5]: size of the BLR 
in units of 10 15 cm; Col. [6]: power injected in the blob calculated in the comoving frame, in units of 10 45 erg s — , Col. [7]: accretion disc luminosity in 
units of 10 45 erg s _1 and (in parenthesis) in units of £-Edd; Col. [8]: magnetic field in Gauss; Col. [9]: bulk Lorentz factor at -Rdiss; Col. [10] viewing angle 
8 V in degrees; Col. [11] and [12]: break and maximum random Lorentz factors of the injected electrons; Col. [13] and [14]: slopes of the injected electron 
distribution [Q(~f)] below and above 7^,; The total X-ray corona luminosity is assumed to be in the range 10-30 per cent of L d . Its spectral shape is assumed 
to be always oc v~ 1 exp(— hu/150 keV). 



Name 


logP r 


logP B 


logP c 


logPp 


0937+5008 


43.49 


42.06 


44.47 


45.05 


1040+23 


43.89 


43.43 


44.08 


44.99 


1055+01 


45.31 


44.69 


44.68 


46.67 


1 106+023 


42.87 


41.76 


44.58 


46.94 


1117+2014 


42.98 


43.66 


42.46 


42.27 


1509+022 


44.42 


42.04 


45.04 


46.03 



Table 7. Logarithm of the jet power in the form of radiation, Poynting flux, 
bulk motion of electrons and protons (assuming one proton per emitting 
electron). Powers are in erg s — . 
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Figure 7. SED of 1055+01 and 1106+023. Lines as in Fig.[6] 



© 0000 RAS, MNRAS 000, 000-000 





© 0000 RAS, MNRAS 000, 000-000 



